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through these geometrical objects. 
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1. Introduction 

Three-level systems are of fundamental importance to many branches of physics. While 
two levels give the simplest model for the dynamics of discrete systems, three levels 
illustrate the role that an intermediate state can play in inducing transitions between 
the other two. Canonical examples of this include applications in quantum optics that 
use three-level atoms to control quantum state evolution [Ij. Such laser control is 
used, for instance, to transfer population between two states using stimulated Raman 
adiabatic passage (STIRAP) [21 [3] and chirped adiabatic passage (CARP) In some 
of these systems, the interaction of the radiation with the atom is represented as a 
time-dependent Hamiltonian inducing an energy separation between the two states 
that varies with time. For a non-zero sweep rate, it can be shown that there is finite 
transition probability between the states [HIE], [7]. The study of Landau-Zener transitions 
in multilevel systems is of interest to understand the interplay between various level 
crossings [8]. Particle physics represents another example where three- level systems 
play a central role as, for example, the oscillations of neutrino flavor eigenstates [9]. 

The general Hamiltonian of a three-level system involves 8 independent operators. 
Such a set can also naturally arise as a subgroup of higher level systems where there 
is some degeneracy involved. Thus, several important two-qubit problems in quantum 
computing and quantum information can be so written in terms of eight operators 
that form a subalgebra of the full fifteen operators that describe two spins. The 
Hamiltonian describing anisotropic spin exchange is an example of one such important 
physical problem. While isotropic spin exchange has been explored to design two-qubit 
gates in quantum computing, anisotropic spin exchange has been studied as a possible 
impediment to two-qubit gate operations [TDl [H]. Such a SU(3) Hamiltonian is given 
by 

H(t) = J{t){a.T + (3{t).{a x r) + ^ .r(t).f), (1) 

when written in terms of a scalar, a vector and a symmetric tensor operator expressed 
in terms of two Pauli spins. Here, f3{t) is the Dzyaloshinksii-Moriya vector [121 [13] and 
r(t) is the (traceless) symmetric interaction term. While the first term is the familiar 
Ising interaction Hamiltonian [T3|, the last two terms are due to spin-orbit coupling. 

Given this wide applicability, a geometrical picture of the dynamics of three-level 
systems can be useful. For a two-level system, the geometry of the evolution operator 
is well known. Any density matrix can be written as p = (J*-^-' -|- n.a')/2, where a are 
the Pauli matrices. Unitary evolution of p is represented as the vector n rotating on 
the surface of a three dimensional unit sphere called the Bloch sphere [15]. This vector, 
along with a phase, accounts for the three parameters describing the time evolution 
operator of a two-level system. The vector n, along with the phase factor, is shown in 
Fig. ([1]). The vector n shown traces out the "base manifold" and together with the 
global phase factor or "fiber" at each point on that manifold is referred to as a "fiber 
bundle" While the density matrix is independent of it, the complete description of 
the system requires this phase as well. The aim of this paper is to provide an analogous 
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geometrical picture for a three-level system with appropriate generalizations of the base 
and fiber. 

Some work already exists regarding the geometry of SU(3). Following Wei and 
Norman [17], Dattoli and Torre have constructed the "Rabi matrix" for a general SU(3) 
unitary evolution in [18]. Mosseri and Dandoloff in [T9| described the generalization 
of the Bloch sphere construction of single qubits to two qubits via the Hopf fibration 
description. This method relies upon the homomorphism between the SU(2) and S0(3) 
groups and likewise between the SU(4) and SO (6) groups. In [20], the authors propose 
a generalized Euler angle parameterization for SU(4). This decomposition is similar to 
the work in [211 ESI El [25l [26l [27] into which fits our treatment of SU(3) in this 
paper. 

Another well known choice of the (A^^ — 1) generators Sj of the SU(N) group 
was studied in [281 [29] . Consider Sj, chosen to be traceless and Hermitian such that 
[sj,Sj] = 2ifijkSk and Tr{sjSj} = 26jk- Here, fijk is the completely antisymmetric 
symbol which for a two-level system is the Levi-Civita symbol eijk, and a repeated index 
is summed over. In this basis, the Hamiltonian is written as H(t) = FjSj. With this 
choice, the Liouville-Von Neumann equation for the density matrix p = I/N + SjSj/2 
becomes Si = fijk^jSk- Note that for the N=2 case, this is the familiar Bloch sphere 
representation. But, for SU(3), this representation differs from the one we present 
in two aspects. Firstly, the "coherence vector", whose elements are real and are 
given by 5*^, experiences rotations in a (A^^ — 1) dimensional space. For instance, 
for SU(3), the coherence vector undergoes rotations in an eight-dimensional space. 
Arbitrary rotations in eight dimensions are characterized by 28 parameters. But since 
a three-level Hamiltonian is only characterized by 8 real quantities, this means that 
the coherence vector is not permitted arbitrary rotations and is instead constrained. 
Secondly, the coherence vector representation does not differentiate between local and 
non-local operations. Our decomposition of the time evolution operator into a diagonal 
and an off-diagonal term in this paper is more suited for this differentiation. Such 
a parameterization of the time evolution operator in terms of local and non-local 
operations can be useful in understanding entanglement. The aim of this paper is 
to discuss the geometry of two-qubit time evolution operators in terms of such a 
decomposition. The authors in ^30] discuss an alternative decomposition of two-qubit 
states in terms of two three- vectors and a 3 x 3 dyadic to discuss entanglement. 

A series of papers presented a systematic approach to studying N-level systems 
using a program of unitary integration [2Tl[22l[23l[2il[3ll[32l[25l[26l[27]. Continuing 
this program, we present a complete analytical solution to the three-level problem 
that generalizes the Bloch sphere approach to three levels. Below, we define the 
fiber bundle via two different decompositions which allows us to extract the geometric 
phases associated with a three-level system (for a discussion on the quantum phases 
of three-level systems, see [331 [Mj)- These fiber bundles are {SU(3)/SU(2)x U(l)} x 
{SU(2)xU(l)} and {SU(4)/[SU(2)xSU(2)]} x {SU(2) xSU(2)}. The structure of this 
paper is as follows: Section [2] outlines the unitary integration program to solve time- 
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Figure 1. Bloch or Poincare sphere representation for SU(2). The base manifold is 
the sphere while the fiber is given by the U(l) phase at each point on that sphere. 
Together, we have the fiber bundle SU(2) ~ 8^x11(1). 



dependent operator equations. Section [3] uses this technique for the solution of a general 
time-dependent SU(3) Hamiltonian completely analytically. Section H] presents the 
geometry of the time evolution operator for SU(3) with some applications. Section 
|5] presents a coordinate description that is useful to define the geometric phase for 
three-level systems, and Section [6] presents the conclusions. The appendix will present 
an alternative analytical solution to the three-level problem by exploiting the natural 
embedding of SU(3) in SU(4). 

2. Unitary Integration 

Many important applications in physics involve time dependence in the Hamiltonian. 
For such systems, the time evolution operator is not given by the simple exponentiation 
of the Hamiltonian [35] . To handle the time evolution for such Hamiltonians iteratively, 
"Unitary Integration" was proposed in [211 [22|, [231 121] ■ Earlier work with this technique 
is presented in [171 EH]. Later, the technique was presented as generalizing the SU(2) 
example to solve iteratively for the time evolution operator \J^^\t) of N-level systems 
|31[ [32] . Consider the N-dimensional Hamiltonian H^^^ given by 



The diagonal blocks are (N— n)- and (n)-dimensional square matrices, respectively, while 
V is an {N — n) x (n)-dimensional matrix. 

The evolution operator U*^^^(t) for such a H^^^ is written as a product of two 
operators U'^^^(t) = U1U2, where 




(2) 




(3) 



Block sphere like construction of SU(3) Hamiltonians using Unitary Integration 



5 



For any N, n is arbitrary with 1 < n < A^, and tilde denotes that the matrices need 
not be unitary. The product of three factors parallels the product of exponentials in 
three Pauli matrices. Equations defining the rectangular matrices z{t) and w^{t) are 
developed and the problem is reduced to the two residual (A^ — n)- and (n) dimensional 
evolution problems sitting as diagonal blocks of 1/2- z{t) and w"''(t) are related to each 
other through the unitarity of \J^'^\t) [311 [32] : 

z = -71 w = -W72, (4) 

with 7i = + z.zt and 72 = i^"^ + z^.z. 

With \J^'^\t) in such a product form, the Schrodinger equation is written as 

tlf^it) = Heff t/2, (5) 

Heff = t/r^H(^)f7i - tU^'U,. 

Here, overdot denotes differentiation with respect to time. Since U2 is block diagonal, 
the off-diagonal blocks of equation (JSj) define the equation satisfied by z 

iz = H(^-")z + V-z(Vtz + H(")). (6) 

Note that the initial condition f/^(0) = implies that f7i(0) = U2{0) = 

and z(0) = 0'^^""^ equation IQ, along with the initial condition can be solved to 
determine z and thereby Ui and Hgff for subsequent solution of equation ([5]) for U2. In 
this manner, the procedure iteratively determines U^^'^{t). 

Before discussing the geometry of the time evolution operators for this unitary case, 
we briefly mention the procedure to deal with non-Hermitian Hamiltonians. For such a 
non-Hermitian Hamiltonian, 

= ( J ] (7) 

where tilde denotes possibly non-Hermitian character, and the off-diagonal components 
V and Y are independent. In this case, equation IQ is replaced by 

iz = H(^-")z + V -z(Y^z + H(")), (8) 
and there is a separate equation governing the evolution of w given by 

iw^ = wt(zYt - h(^-")) + (hW + Ytz)wt + Y^ (9) 
The diagonal terms of the time-evolution operators are governed by 

lU2{t) =i ~,, ,]U2. (10) 

Returning to the case where the Hamiltonian is Hermitian, it is convenient to 
render the two matrices Ui and U2 themselves unitary [311 [32]. For this purpose, a 
"gauge factor" b is chosen such that the unitary counterparts of Ui and U2 are defined 
via Ui = Uib and U2 = h~^Ui. Since UllJi = diag{'y[ ^^72), this would imply that b is 
the "Hermitian square-root" of diag{'y[~^\'y2)- This "Hermitian square-root" is defined 
by the relation [b^~^^yb^~^^ = diag{'~f[ ^^72)- Inspection of the power series expansion 
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of 7i^^'' = (i + z.z^)*^^^) and 72^^^ = (I + z^z)'-^^) show that since each term in the 
expansion is Hermitian, matrices 7^ ^ and 72 ^ are Hermitian and have non-negative 
eigenvalues. Because of this, it is sufficient to define h as the inverse square root via 

Furthermore, H^ff in equation ([5]) is Hermitian for the unitary counterpart Ui. The 

upper diagonal block of this Hermitian Hamiltonian accompanying the decomposition 

U = U1U2 is given by 
_ 1 

^[%^,7f ] + I (7^(5^"-'^) - zVt)7f + H.c) , (11) 

where [,] represents the commutator and H.c. stands for the Hermitian conjugate. The 
lower diagonal block is similarly given by 

_ 1 

^[^,7|] + I [if^iHi---^ + ,tv)^| + H.c) . (12) 

For = 3, n = 1, these diagonal blocks define an SU(2)- and a U(l) Hamiltonian 
and z is a pair of complex numbers. The SU(2) Hamiltonian is in turn rendered in 
terms of its fiber bundle in Fig. ([T]) and the U(l) Hamiltonian corresponds to a phase. 
Together, they describe a four- dimensional fiber for SU(3) over the base manifold, also 
four dimensional, of z. 

Alternatively, = 3 SU(3) problems may be conveniently seen as a part of = 4 
SU(4) problems, making contact with two qubit systems that are extensively studied. 
In this case, for = 4, = 2, these diagonal blocks define two SU(2) Hamiltonians and 
z is a 2 X 2 matrix representable in terms of Pauli spinors. Generally, it is 8-dimensional 
while the fiber has seven dimensions (two SU(2) and a mutual phase) but for the SU(3) 
subgroup of SU(4),both the base and manifold again reduce to four dimensions each. 
With z a pair of complex numbers, the non-trivial part of geometrizing SU(3) is thereby 
reduced to describing this four- dimensional manifold. Exploring this for the A^ = 3, 
n = 1 decomposition will be the content of the next section whereas the Appendix gives 
the alternative SU(4) rendering. 

3. Geometry of general SU(3) time evolution operator 

A general time-dependent three-level Hamiltonian may be written in terms of eight 
linearly independent operators of a three-level system. Such a Hamiltonian can also 
be written in terms of a subgroup of 15 operators of a four-level system. Before the 
time evolution operator is presented in the SU(3) basis in terms of a A^ = 3, n = 1 
decomposition, we will note that it can be rendered in a few alternative ways. 

First, a general time- dependent four-level Hamiltonian may be written as H{t) = 
X^jCjOj. Here Cj are time-dependent and Oj are the unit matrix and 15 linearly 
independent operators of a 4-level system that may be chosen in a variety of matrix 
representations. One choice used in particle physics are the so called Greiner matrices 
[361 [2n [22I [231 [2i] . Another choice consists of using ct, r, <t ® r and the 4x4 unit 
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matrix. Such a choice was discussed in [23 [2S] and will be used throughout this paper. 
As it stands, the above Hamiltonian describes a general four-level atom with 4 energies 
and 6 complex couplings. Note that only the three differences in energies are important. 
Restricting the 15 coefficients q to a smaller number allows this Hamiltonian to describe 
various physical Hamiltonians, forming different subalgebras of the su(4) algebra |25j . 
For example, if two of the six complex couplings are zero (levels 1 and 4 and levels 2 
and 3 of a four- level atom not coupled), then the Hamiltonian may be recast such that 
the operators involved belong to an so(5) subalgebra [25]. On the other hand, if levels 
2 and 3 are degenerate and level 4 is uncoupled from the rest, then the problem may be 
recast in terms of only eight operators belonging to the su(3) subalgebra of su(4). This 
is illustrated in Fig. [2] and is one of the systems of interest in this paper. 




Figure 2. Levels |2) and |3) couple equally to |1) and to |4), which are themselves 
coupled. The three complex coupling matrix elements and two energy positions define 
such an SU(3) system. 

Alternatively, after one arrives at the linear equation for the = 4, n = 2 
decomposition, one can represent the resulting vector in terms of six homogeneous 
coordinates. This is the so-called "Pliicker coordinate" representation for the SU(3) 
Hamiltonian. These coordinates as well as the alternative derivation are presented in 
the appendix. The N = 3, n = 1 decomposition will be the content of the rest of this 
section. 

Consider the Hamiltonian in the basis of the Gell-Mann lambda matrices [37] 
H{t) = J2i The A^ = 3, n = 1 decomposition consists of writing the time evolution 
operator in terms of a product of two matrices U = U1U2 where Ui is composed of a 
(2xl)-dimensional z, as explained in Sec. II. The equation that governs the evolution 
of z, equation ([6]), can be written in this 

= -iVf, - iF^^z^ + iV*z^Zf,; = 1, 2. (13) 

Here, the symbols used in defining z are defined as V = (04 — ia^, ag — ia-i), and 

_ ( ^3 + V^as ai - \ 
\ ai + ia2 —as + v^ag / 
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Using the transformation equations 2 = — 2:1,2 (-De**^)"^, = (De*'^)"^ and + 
|m2p + |m3p = 1 leads to the evolution equation for m = {rriir, rn2r, m^r, rriu, m2i, m^i)^: 
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which describes the rotation of a unit vector in a six dimensional space of the real 
and imaginary parts of fh defined by = m^r + ^^^i- In the above equations, 
D = {1 + l^ip + 1-22^)^^^ and i(j) = i{V*Zi, + Vuzl). The phase is real and determined 
only up to a constant factor. Since the real and imaginary parts of mg are not 
independently defined, the geometrical description of the base manifold for the = 3, 
n = 1 decomposition may be thought of as a point on the surface of a constrained 
six-dimensional unit sphere. 

The two constraints, namely |mip-|- |m2p + |m3p = 1 and the "phase arbitrariness" 
of 0, reduce the 6-dimensional manifold of the three-dimensional complex vector fh 
to a four- dimensional manifold in agreement with there being only four independent 
parameters in z.The first condition defines the base as a vector on an sphere while 
the phase arbitrariness serves as an additional constraint. The fiber, on the other hand, 
is an SU(2) block, evolving as a vector on S'^ Poincare-like sphere with a phase at each 
point, and a U(l) block that amounts to an extra phase. This is presented schematically 
in Fig.Q, as the product of three matrices of the evolution operator. 
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Figure 3. The base and fiber for tlie SU(3) group. Tlic first two factors give tlie 
base manifold, an spliere witfi a phase arbitrariness defined in the text. The fiber, 
described by the third matrix, is composed of a Bloch sphere and a phase associated 
with each of its points, and the second an extra phase represented by a vertical line. 

The alternative = 4, n = 2 decomposition in the appendix yields the equation of 
motion for = —z^/De^'^ in equation ( lA.llI) . Following equation ( ITTi) and equation 
( IT2|) . we see that for this case, the two remaining blocks of the time evolution operator, 
namely t/*^^"^-* and U^'^\ can be transformed into unitary matrices for SU(2). The fiber 



Block sphere like construction of SU(3) Hamiltonians using Unitary Integration 



9 



evolves as vectors on two identical 5*^ Bloch-like spheres with a mutual phase, whose 
evolution is coupled to the base that evolves as a vector on an 5*^ sphere. This is 
illustrated in Fig. (jl]). Either decomposition can be used to study various physical 




Figure 4. The base and fiber for tlie SU(3) group via the iV = 4, n = 2 decomposition. 
The base again is given by an sphere as in Fig. ([3]). The fiber is composed of two 
identical SU(2) Bloch spheres plus phase, and an extra mutual phase between them. 
The four parameters each of base and fiber again account for all eight parameters of 
the SU(3). 



processes as will be discussed in the next section. 



4. Applications 

It is often desirable to control the time evolution of quantum states to manipulate 
an input state into a desirable output state. In [38l [39], the authors considered a 
Hamiltonian of the form Hq — fJ'S{t), where Hq is a free-field Hamiltonian and ^£{t) 
is a control field. To illustrate the "Hamiltonian encoding" scheme to control quantum 
systems, the authors considered a three-level system and studied stimulated Raman 
adiabatic passage (STIRAP), an atomic coherence effect that employs interference 
between quantum states to transfer population completely from a given initial state to a 
specific final state. This is done through a "counterintuitive" pulse sequence. Consider 
the Hamiltonian 

/ Giit) \ 
H(t)= Gi{t) 2A G2{t) . (15) 
V G2{t) ] 

Here G'i,2(i) = 2.5exp[— (t — ^1,2)^/^^] and A = 0.1. The initial population is in the 
upper state. For ti = r, t2 = and r = 3, it is seen that the two empty states are 
coupled first via 6*2 (^) and then the levels |1) and |2) are coupled through Gi. The 
dynamics of the populations reveal complete population transfer. A complete solution 
as per Section [3] was constructed for this model and the results are presented in Fig. [5] 
in total agreement with the results of [38] • 
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Figure 5. Population Pij = plotted as a function of time. The initial 

population in state |1) is completely transferred to |3). Both the unitary integration 
solution and the direct numerical solution [38] are plotted and they coincide at all 
times. 




Quantum control can also be achieved by understanding the nature of tunneling. 
The famous Landau-Zener formula O El [7] predicts the transition probability of the 
ground state of a two-level system when the energy levels adiabatically undergo a 
crossing. The study of level crossings has since been extended to multi-level systems. 
For example, in [ID], the authors considered a three- level atom to study population 
trapping by manipulating the phase acquired as a three-level system evolves under the 
influence of frequency modulated fields [H]. Such a frequency modulated field is given 
by 

E{t) = Eie-^['^i*+'^iWl + E2e-^['^2*+'^2Wl + c.c. (16) 
<^i{t) = Mi sin flit. (17) 

Here, c.c. stands for complex conjugation. The phase (fi{t) in the exponent can be 
written in terms of Bessel functions as [12] 

oo 

gAf,sinC,t_ ^ Jfc(Mj)e*^^^*. (18) 
fc=— oo 

For large values of Qj, the leading contribution for slow time scales would come from 
Jo{Mj). Hence, for large flj, the interaction Hamiltonian can be written as 

H,„i(t) = -d.(Ei Jo(Mi) + E2MM2). (19) 

Hence, for values of Mi^2 that are zeros of the zeroth-order Bessel functions, the 
interaction Hamiltonian is zero and population trapping is observed. Under this 
assumption, consider the full Hamiltonian under the rotating wave approximation, 

/ E,{t) Gi(t) \ 
H(t)= Glit) G2{t) . 

[ cm E,{t) ) 

Here, Ei{t) = Ai - MiQi cos(fiit + 6) and E^lt) = -A2 + M2Q2 cos(fi2t). Results are 
presented in Fig. [6], and for the parameter values ^1,2 = 1, Ai = — A2 = 10, 6* = and 
Gi,2 = 6, demonstrate the phenomenon of population localization discussed in [iO] . 
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Figure 6. (a) For Mi^2 = 7 and the other parameter values given in the text, there is no 
population trapping observed, (b) The energy landscape for Mi 2 — 30.6346 showing 
energy level crossing, (c) Population trapping is observed with Afi,2 = 30.6346 which 
corresponds to the tenth zero of the zeroth-order Bessel function. Note that the thick 
line is Pu and the thin line corresponds to P12. The results agree completely with 



As a final illustration of the unitary integration technique applied to three-level 
systems, let us consider the example discussed in |13]. Here, a three- level system is 
subject to strong fields and the correlation between the scattered light spectrum and 
the atom dynamics is discussed. The authors consider the Hamiltonian 



H(t) 



Gi{t) \ 
G2{t) 




(20) 






\G\{t) cm 

Here, Gi,2(^) = — Vi,2e~*''*. The time evolution of the states calculated as per our 
procedure in Section [3] is plotted in Fig. [7| for different values of the parameters. All of 
these results agree with those given in [l3j. Further features of the base and fiber will 
be presented at the end of the next section. 
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Figure 7. (a) Populations Pij = for S = 5, Vi = 2 and V2 = 1. Pii is given 

by the solid line and P12 is given by the thin line, (b) Same as (a), for S = 12. Note 
that Pi3 oscillates close to zero at all times, (c) Pij for 6 — 12, Vi — 1 and V2 — 2. 



5. Geometric phase for SU(3) group 

Many physical systems give rise to a measurable phase that does not depend directly 
on the dynamical equations that govern the evolution of the system, but depends only 
on the geometry of the path traversed by vectors characterizing the state of the system. 
This geometric phase is denoted by ■jg and is given by the integral [H], 

7, = 1 t/R .(n(R(t))|«VR|n(R(t))), (21) 

where the state evolution is governed by a set of internal coordinates that parameterize 
the Hamiltonian R(t), and Vr is the gradient in the space of these internal coordinates. 
This phase has been generalized to non-cyclic non-adiabatic evolution of quantum 
systems [15| B6| HTj. The purpose of this section is to present this phase in terms 
of coordinates on the Bloch sphere for two-level systems and extend it to three-level 
systems. 

In two-level systems, the time evolution operator is described by three parameters 
as described in Section [H Two of these parameters describe a point on the Bloch 
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sphere. Traversing closed loops on this Bloch sphere returns the quantum system to its 
initial state as described by the two parameters on the Bloch sphere but not the third 
parameter of an overall phase. Hence, general closed loops on the Bloch sphere do not 
correspond to closed loops in the space of the full unitary operator. This discrepancy in 
the phase between the initial and final state corresponds to the geometric phase given 
above and amounts to changes along the fiber at each point on the sphere. To formalize 
this, consider Ui, given by equation ([3]) as unitarized through the matrix b in Section [2l 
which for = 2, n = 1 takes the form 



By identifying cos | = (1 + |zp)~^ and sin |e~*'' = — z(l + |zp)^5j we get the usual 
description of the base manifold in terms of the angles < 6 < n and < e < 27r that 
are associated with the Bloch sphere, namely. 

In terms of the parameters 6 and e, the Hamiltonian H{t) = —a. a is given by 

H{t)=( -^i-^f^y (24) 

^ ' V - sin 9e'' cos 9 J ^ ^ 

equation governing the evolution of the fiber U2 has two terms. The first term is 
evaluated as 

UlH{t)U,= (^-^^ (25) 

which corresponds to the eigenvalues of the Hamiltonian. To evaluate the second term, 
consider the case whereby the vector on the Bloch sphere traverses a closed path defined 
by a constant 6. The second term is then given by 

rrl dU, _ ( -sin^f -|sin^e-M 

Integrating e from to 2tt yields 

f'\.rri dU, /7r(l-cosg) \ 

h ^d{-it) \ -7r(l -COS0) ) ' 

which is the correct formula for the geometric phase of a two-level system |44j . 

To extend this analysis to three-level systems, we consider the = 3, n = 1 
decomposition. The matrix Ui = Ui.b is now given by 



/(2) 



rZZ 




U, = I ^(^+^) ^ I , (27) 



where z is a complex column vector {zi, 2:2)^ and D = yl -\- |zp. Care has to be taken 
in assigning angles to elements of this matrix such that the transformation satisfies 
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two conditions: the Ui matrix should not depend on and the transformation must 
be commensurate with the definition of fh. To this effect, we transform z into polar 
coordinates: Zi = — tan y cos %-e*'^% Z2 = — tan ^ sin ^e**^^. These transformation 
equations imply that D = 



1+ z 



sin ^ sin ^e^^"^- 



mi 



( 



\ 



1 

sm 



and m3 = cos ye The Ui matrix is given by 
2^ -sin^ ^sin^2e*('^"'2) 

i{e-l-f-2) 



sin Y cos yc**^^^" 



1112 



- 2 sin^ ^ cos^ f 
2 ^ sin e-).e 



sin f cos f e-"i 



sin ^ sin 6'2e* 
1-2 sin^ ^ sin^ f 



sin Y sin ye" 



sin Y cos ye' 
sin Y sin yC*' 



cosf 



(2^ 



In the above equation, the range on the angles < < vr and < e-j < 27r 
are chosen so that the absolute value of each element of the time-evolution operator is 
positive [IE]. Hence Ui can be represented as two vectors on a sphere, at angles (6'i, ei) 
and (6*2, £2) respectively. This is represented in Fig. ([8]). Since the columns of a unitary 




Figure 8. The base manifold Ui is characterized by two sets of angles < 6',; < tt, 
< e.; < 27r which can be represented as two vectors with polar angles (6*1 , ei ) and 

(02, £2). 



operator correspond to normalized eigenvectors, we can consider the last column of 
the matrix above, = (— sin y cos ^6*"^, — sin y sin yc^^^, cos y)"^, and evaluate the 
so-called connection 1-form given by [19] 

A=-i{^\d\ij). (29) 

The Abelian geometric phase, given by 7^ = / ^ is evaluated to be 

If 9 

7g = -- y sin^ — {{dei + dt2) + cos02(c?ei - dt2)). (30) 

If the various angles are relabelled ei — *• —7 — a, €2 — ^ —7 + a, 61 ^ 29 and 6*2 2j3, 
the formula above agrees with [50] and [H] . The time-evolution operator above can now 
be used as in the case of SU(2) to evaluate the dynamic contribution / UlH{t)Ui and 
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the geometric contribution to the time evolution operator which is given by —i J UldUi, 
where dUi = ^dOi + ^de^ z = 1, 2. 

This description of the base manifold in terms of {6i, €{) can now be used to describe 
the dynamics of various physical processes. Fig. ([9]) represents the base manifold 
corresponding to the results in Fig. ([71). (6'i,ei) depend on all the parameters that 
define the system while (6*2,62) depend only on the ratio V1/V2. Also note that the 
maximum value of €2, corresponding to the maximum latitude traversed by the black 
curve, is inversely proportional to 6. Such observations can be used to control the 
dynamics of this system. 




Figure 9. The base manifold corresponding to tlie results in Fig. ((7|) for the three- 
level system of [43^. For the first column, Vi = 1, V2 = 2. The second column 
corresponds to T^i = 2, V2 = 2 and the third to Vi = 2, V2 = 1- The rows correspond 
to 5 = 1, S — 5 and 6 — 50. The thin black curve describes {9i, ei) and the thick red 
curve the set (6*2, £2)- 

6. Conclusions 

The ability to decouple the time dependence of operator equations from the non- 
commuting nature of the operators is the central feature of unitary integration and 
also characterizes the Bloch sphere representation for the evolution of a single spin. By 
doing so, the quantum mechanical evolution is rendered a "classical" picture of a rotating 
unit vector. For a two-level atom, the Bloch sphere representation along with a phase 
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completely determines the time evolution operator. In this paper, we have extended 
this program to deal with the time evolution operator belonging to the SU(3) group. 
This complements the work in [3l] for SU(4) Hamiltonians of two qubit systems. We 
have also extended the analysis of geometric phase to three-level systems by providing 
an explicit coordinate representation for the SU(3) time evolution operator. 



Appendix A. Alternative derivations for a general SU(3) Hamiltonian. 

Consider a three-level Hamiltonian written in terms of the Gell-Mann matrices [37] as 
H(t) = J^i^iaiXi. To exploit the fact that this Hamiltonian is a subgroup of four-level 
problems, it is represented in terms of the O matrices [25] as 

2-^02 + (ag - ^)03 + (2a3 + 2-^)04 + 0405 + a^Oe + ^a^O-j + 2a508 + 

aiOg + asOio + 2aiOii + 2a20i2 + 2a60i3 + 2a60i4 - 2ajOi^ + 2a^Oi&. (A.l) 

This embeds the Hamiltonian H{t) = J^i^iXi as a 4x4 matrix with zeros along the 
last row and column. In such a representation, the various entries of the Hamiltonian 
equation ([2]) are given by 

H'^^'^^ = ^agl^^^ + aiCTi + 020-2 + 030-3, (A.2) 
v3 

V = ^(04 - ^5)1(2) + l(a6 - ^7)0-1 (A.4) 

Writing z in the standard Clifford basis as z = ^z^I^'^^ — ^Y^i^iCTi, it follows from 
equation ([6]) that zi = iz2 and z^ = iz4^ and the equation reduces precisely to equation 
(IT^ . The geometry described in Section [3] can thus be derived from either of these 
decompositions of the time evolution operator. 

The SU(3) subgroup in equation (lA.ip is one among many SU(3) subgroups 
embedded in SU(4). Another choice corresponds to the Dzyaloshinskii-Moriya 
interaction Hamiltonian [121 113] and is also of interest because the 4x4 matrices now 
do not have a trivial row and column of zeros. In the two-spin basis, this Hamiltonian 
is given by 

Hit) = qO, = ai(02 + O3) + 2a2(Oi5 + Oie) + 2a3(Oi4 - O13) + 204(07 + On) 

i 

+05(05 + Oio) + 05(05 + O9) + 207(03 + O12) + ^(204 - Oi3 - O14). (A.5) 

The correspondence between the coefficients in terms of O and in terms of the A matrices 
is : ci = 0, C2 = oi, C3 = oi, C4 = 4o8/ a/3, C5 = Og, Cq = 05, C7 = 204, Cg = 207, Cg = oe, 
Cio = 05, cii = 204, C12 = 207, Ci3 = -203 - 2as/V5, Ci4 = 203 - 2o8/\/3, C15 = 2o2 
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and Ci6 = 2a2. Relabeling of the states 1— s>2, 2— s>3, 3-h>4 and 4 — > 1 expresses the 
Hamiltonian as 

-ff^"^"^^ = ^agl^^^ - a^cTi - a2(T2 - aia^, (A. 6) 

v3 

H^'^ = -^«8l(^) - ^asa„ (A.7) 
V = i(a6 - ^7)1(2) + l(a6 - ^7)0-1 (A.8) 

-^(os + «a4)cr2 - -(04 - ^asjcTs- 

If z is written in terms of the standard Clifford basis (I, —ia) as z = ^Z4l^'^^ — | I]^=i -^^i'^'i 
, it follows from equation ([6]) that zi = iz^ and ^2 = iz^- This is consistent with the 
parameter count that since the inhomogeneity V has only two free complex parameters 
(namely Vi = — ia-j and V2 = (24 — ia^)-, the complex z matrix should be composed 
only of two independent complex parameters, Zi and Z2. With the above analysis, 
equation ([6]) becomes for the pair of complex numbers 

]^z^ = ^Xf, - iF^^Zy + 2GyZuZ^] fi,u = 1, 2. (A. 9) 

Here X = {Vi/2, -iV2/2), G = {2V*,2iV2*) and 

■p _ ( ^^3 — V^^Cts 0-1 + ^0-2 \ 

\ —ai + ia2 —ia-i — \/^ia^ ) 
Paralleling the technique employed to solve an SO (5) Hamiltonian in [311 132] , 
we transform z into a complex vector m: = and ma = ^ such that 

|mip + |m2p + |m3p = l,with D = {1 + ^{\zi\^ + |z2n)"'"^^- This leads to the new 
set of evolution equations 

ias — ^f^ia^ oi + ia2 — og + ^a^ \ 
—ai + ia2 —ias — y/Sias as + ia4 m. (A. 10) 

ag + ia^ —as + ^04 J 

This can be written as an equation describing the rotation of the real and imaginary 

components of the vector fh = (mi^, m2r, m^r, f^u, fn2i, fT^u)^ 



m 



1 



m 



/ 





ai 


-ae 


—as + VSas 


— a2 


— a7 


\ 




— ai 





as 


—a2 


03 + V3a8 


— a4 








-as 





—aj 


— a4 











a2 


07 





ax 


-ag 






a2 


-as - y^as 


a4 


-ai 





as 




V 


aj 


04 





ag 


-as 





/ 



rh. (A.ll) 



Here, the coefficients Cj are written in terms of the coefficients a^, whose correspondence 
was given earlier in this section. Also note that = m^^ + D = (1 + |-2;2p)^ 
and = {V*Zi, + V^z*). Simplifying this leads to the equation i(p = — 2(X^z* — X*^z^) 
for the evolution of which is clearly real but determined only to within a constant. A 
little algebra yields for the effective Hamiltonian given by equation (ITT]) . 
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and for the effective Hamiltonian given by equation (1121) . the expression H^'^^ + (z^V + 
Vtz)/2. 

Another representation of the SU(3) subgroup of SU(4) Hamihonians is given by the 
so called "Pliicker coordinate" representation of the SU(4) group discussed in [311 [32] . 
For an arbitrary SU(4) matrix, the Pliicker coordinates are defined as a set of six 
parameters (P12, P13, P14, -P23, -^24, -P34) such that P12-P34 - -P13-P24 + -P14-P23 = and 



E \P^, 



1. They can be written in terms of the unit vector m and are given by 



( P12 






/ 


irriQ — ms 


\ 


Pl3 








irrii + m2 




PlA 




1 




—irris + m4 




P23 




~ 2 




—im^ — m4 




P24 








—irrii + m2 




[ P34 


J 




I 


irriQ + ms 


/ 



(A.12) 



The linear equation of motion for m translates into an evolution equation for P 
(P12, -Pi3, -P14, -P23, -P24, -P34) of the form iP = HpP. Here, Hp is given by 



H, 



H 

Vl 



PI 




(A. 13) 



where 



H 



H 



PI 



P2 



V, 



/ 2a8/y3 
y 064- - ia75- 
-a8/V3 

064- + ^075- 

^ —(364+ — ^0-75+ 
^32- 

V 



064 



V 



+ ^075 

-08/^3 

ai 

-064- + «a75 
(364+ + ^^75+ 



-^32- 



064- + ^a^5- 
og/VS 
-ai 



-(364- — ?Ct75- 
-(364- — ■^075- 

-208/73 







— 064^ 
064+ 



- ^075 
2075 



+ / 



In the above equation, Ojj± denotes Oi ± Oj. 
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1. Introduction 

Three-level systems are of fundamental importance to many branches of physics. While 
two levels give the simplest model for the dynamics of discrete systems, three levels 
illustrate the role that an intermediate state can play in inducing transitions between 
the other two. Canonical examples of this include applications in quantum optics that 
use three-level atoms to control quantum state evolution [Ij. Such laser control is 
used, for instance, to transfer population between two states using stimulated Raman 
adiabatic passage (STIRAP) [21 [3] and chirped adiabatic passage (CARP) In some 
of these systems, the interaction of the radiation with the atom is represented as a 
time-dependent Hamiltonian inducing an energy separation between the two states 
that varies with time. For a non-zero sweep rate, it can be shown that there is finite 
transition probability between the states [HIE], [7]. The study of Landau-Zener transitions 
in multilevel systems is of interest to understand the interplay between various level 
crossings [8]. Particle physics represents another example where three- level systems 
play a central role as, for example, the oscillations of neutrino flavor eigenstates [9]. 

The general Hamiltonian of a three-level system involves 8 independent operators. 
Such a set can also naturally arise as a subgroup of higher level systems where there 
is some degeneracy involved. Thus, several important two-qubit problems in quantum 
computing and quantum information can be so written in terms of eight operators 
that form a subalgebra of the full fifteen operators that describe two spins. The 
Hamiltonian describing anisotropic spin exchange is an example of one such important 
physical problem. While isotropic spin exchange has been explored to design two-qubit 
gates in quantum computing, anisotropic spin exchange has been studied as a possible 
impediment to two-qubit gate operations [TDl [H]. Such a SU(3) Hamiltonian is given 
by 

H(t) = J{t){a.T + (3{t).{a x r) + ^ .r(t).f), (1) 

when written in terms of a scalar, a vector and a symmetric tensor operator expressed 
in terms of two Pauli spins. Here, f3{t) is the Dzyaloshinksii-Moriya vector [121 [13] and 
r(t) is the (traceless) symmetric interaction term. While the first term is the familiar 
Ising interaction Hamiltonian [T3|, the last two terms are due to spin-orbit coupling. 

Given this wide applicability, a geometrical picture of the dynamics of three-level 
systems can be useful. For a two-level system, the geometry of the evolution operator 
is well known. Any density matrix can be written as p = (J*-^-' -|- n.a')/2, where a are 
the Pauli matrices. Unitary evolution of p is represented as the vector n rotating on 
the surface of a three dimensional unit sphere called the Bloch sphere [15]. This vector, 
along with a phase, accounts for the three parameters describing the time evolution 
operator of a two-level system. The vector n, along with the phase factor, is shown in 
Fig. ([1]). The vector n shown traces out the "base manifold" and together with the 
global phase factor or "fiber" at each point on that manifold is referred to as a "fiber 
bundle" While the density matrix is independent of it, the complete description of 
the system requires this phase as well. The aim of this paper is to provide an analogous 
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geometrical picture for a three-level system with appropriate generalizations of the base 
and fiber. 

Some work already exists regarding the geometry of SU(3). Following Wei and 
Norman [17], Dattoli and Torre have constructed the "Rabi matrix" for a general SU(3) 
unitary evolution in [18]. Mosseri and Dandoloff in [T9| described the generalization 
of the Bloch sphere construction of single qubits to two qubits via the Hopf fibration 
description. This method relies upon the homomorphism between the SU(2) and S0(3) 
groups and likewise between the SU(4) and SO (6) groups. In [20], the authors propose 
a generalized Euler angle parameterization for SU(4). This decomposition is similar to 
the work in [211 ESI El [25l [26l [27] into which fits our treatment of SU(3) in this 
paper. 

Another well known choice of the (A^^ — 1) generators Sj of the SU(N) group 
was studied in [281 [29] . Consider Sj, chosen to be traceless and Hermitian such that 
[sj,Sj] = 2ifijkSk and Tr{sjSj} = 26jk- Here, fijk is the completely antisymmetric 
symbol which for a two-level system is the Levi-Civita symbol eijk, and a repeated index 
is summed over. In this basis, the Hamiltonian is written as H(t) = FjSj. With this 
choice, the Liouville-Von Neumann equation for the density matrix p = I/N + SjSj/2 
becomes Si = fijk^jSk- Note that for the N=2 case, this is the familiar Bloch sphere 
representation. But, for SU(3), this representation differs from the one we present 
in two aspects. Firstly, the "coherence vector", whose elements are real and are 
given by 5*^, experiences rotations in a (A^^ — 1) dimensional space. For instance, 
for SU(3), the coherence vector undergoes rotations in an eight-dimensional space. 
Arbitrary rotations in eight dimensions are characterized by 28 parameters. But since 
a three-level Hamiltonian is only characterized by 8 real quantities, this means that 
the coherence vector is not permitted arbitrary rotations and is instead constrained. 
Secondly, the coherence vector representation does not differentiate between local and 
non-local operations. Our decomposition of the time evolution operator into a diagonal 
and an off-diagonal term in this paper is more suited for this differentiation. Such 
a parameterization of the time evolution operator in terms of local and non-local 
operations can be useful in understanding entanglement. The aim of this paper is 
to discuss the geometry of two-qubit time evolution operators in terms of such a 
decomposition. The authors in ^30] discuss an alternative decomposition of two-qubit 
states in terms of two three- vectors and a 3 x 3 dyadic to discuss entanglement. 

A series of papers presented a systematic approach to studying N-level systems 
using a program of unitary integration [2Tl[22l[23l[2il[3ll[32l[25l[26l[27]. Continuing 
this program, we present a complete analytical solution to the three-level problem 
that generalizes the Bloch sphere approach to three levels. Below, we define the 
fiber bundle via two different decompositions which allows us to extract the geometric 
phases associated with a three-level system (for a discussion on the quantum phases 
of three-level systems, see [331 [Mj)- These fiber bundles are {SU(3)/SU(2)x U(l)} x 
{SU(2)xU(l)} and {SU(4)/[SU(2)xSU(2)]} x {SU(2) xSU(2)}. The structure of this 
paper is as follows: Section [2] outlines the unitary integration program to solve time- 
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Figure 1. Bloch or Poincare sphere representation for SU(2). The base manifold is 
the sphere while the fiber is given by the U(l) phase at each point on that sphere. 
Together, we have the fiber bundle SU(2) ~ 8^x11(1). 



dependent operator equations. Section [3] uses this technique for the solution of a general 
time-dependent SU(3) Hamiltonian completely analytically. Section H] presents the 
geometry of the time evolution operator for SU(3) with some applications. Section 
|5] presents a coordinate description that is useful to define the geometric phase for 
three-level systems, and Section [6] presents the conclusions. The appendix will present 
an alternative analytical solution to the three-level problem by exploiting the natural 
embedding of SU(3) in SU(4). 

2. Unitary Integration 

Many important applications in physics involve time dependence in the Hamiltonian. 
For such systems, the time evolution operator is not given by the simple exponentiation 
of the Hamiltonian [35] . To handle the time evolution for such Hamiltonians iteratively, 
"Unitary Integration" was proposed in [211 [22|, [231 121] ■ Earlier work with this technique 
is presented in [171 EH]. Later, the technique was presented as generalizing the SU(2) 
example to solve iteratively for the time evolution operator \J^^\t) of N-level systems 
|31[ [32] . Consider the N-dimensional Hamiltonian H^^^ given by 



The diagonal blocks are (N— n)- and (n)-dimensional square matrices, respectively, while 
V is an {N — n) x (n)-dimensional matrix. 

The evolution operator U*^^^(t) for such a H^^^ is written as a product of two 
operators U'^^^(t) = U1U2, where 




(2) 




(3) 
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For any N, n is arbitrary with 1 < n < A^, and tilde denotes that the matrices need 
not be unitary. The product of three factors parallels the product of exponentials in 
three Pauli matrices. Equations defining the rectangular matrices z{t) and w^{t) are 
developed and the problem is reduced to the two residual (A^ — n)- and (n) dimensional 
evolution problems sitting as diagonal blocks of 1/2- z{t) and w"''(t) are related to each 
other through the unitarity of \J^'^\t) [311 [32] : 

z = -71 w = -W72, (4) 

with 7i = + z.zt and 72 = i^"^ + z^.z. 

With \J^'^\t) in such a product form, the Schrodinger equation is written as 

tlf^it) = Heff t/2, (5) 

Heff = t/r^H(^)f7i - tU^'U,. 

Here, overdot denotes differentiation with respect to time. Since U2 is block diagonal, 
the off-diagonal blocks of equation (JSj) define the equation satisfied by z 

iz = H(^-")z + V-z(Vtz + H(")). (6) 

Note that the initial condition f/^(0) = implies that f7i(0) = U2{0) = 

and z(0) = 0'^^""^ equation IQ, along with the initial condition can be solved to 
determine z and thereby Ui and Hgff for subsequent solution of equation ([5]) for U2. In 
this manner, the procedure iteratively determines U^^'^{t). 

Before discussing the geometry of the time evolution operators for this unitary case, 
we briefly mention the procedure to deal with non-Hermitian Hamiltonians. For such a 
non-Hermitian Hamiltonian, 

= ( J ] (7) 

where tilde denotes possibly non-Hermitian character, and the off-diagonal components 
V and Y are independent. In this case, equation IQ is replaced by 

iz = H(^-")z + V -z(Y^z + H(")), (8) 
and there is a separate equation governing the evolution of w given by 

iw^ = wt(zYt - h(^-")) + (hW + Ytz)wt + Y^ (9) 
The diagonal terms of the time-evolution operators are governed by 

lU2{t) =i ~,, ,]U2. (10) 

Returning to the case where the Hamiltonian is Hermitian, it is convenient to 
render the two matrices Ui and U2 themselves unitary [311 [32]. For this purpose, a 
"gauge factor" b is chosen such that the unitary counterparts of Ui and U2 are defined 
via Ui = Uib and U2 = h~^Ui. Since UllJi = diag{'y[ ^^72), this would imply that b is 
the "Hermitian square-root" of diag{'y[~^\'y2)- This "Hermitian square-root" is defined 
by the relation [b^~^^yb^~^^ = diag{'~f[ ^^72)- Inspection of the power series expansion 
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of 7i^^'' = (i + z.z^)*^^^) and 72^^^ = (I + z^z)'-^^) show that since each term in the 
expansion is Hermitian, matrices 7^ ^ and 72 ^ are Hermitian and have non-negative 
eigenvalues. Because of this, it is sufficient to define h as the inverse square root via 

Furthermore, H^ff in equation ([5]) is Hermitian for the unitary counterpart Ui. The 

upper diagonal block of this Hermitian Hamiltonian accompanying the decomposition 

U = U1U2 is given by 
_ 1 

^[%^,7f ] + I (7^(5^"-'^) - zVt)7f + H.c) , (11) 

where [,] represents the commutator and H.c. stands for the Hermitian conjugate. The 
lower diagonal block is similarly given by 

_ 1 

^[^,7|] + I [if^iHi---^ + ,tv)^| + H.c) . (12) 

For = 3, n = 1, these diagonal blocks define an SU(2)- and a U(l) Hamiltonian 
and z is a pair of complex numbers. The SU(2) Hamiltonian is in turn rendered in 
terms of its fiber bundle in Fig. ([T]) and the U(l) Hamiltonian corresponds to a phase. 
Together, they describe a four- dimensional fiber for SU(3) over the base manifold, also 
four dimensional, of z. 

Alternatively, = 3 SU(3) problems may be conveniently seen as a part of = 4 
SU(4) problems, making contact with two qubit systems that are extensively studied. 
In this case, for = 4, = 2, these diagonal blocks define two SU(2) Hamiltonians and 
z is a 2 X 2 matrix representable in terms of Pauli spinors. Generally, it is 8-dimensional 
while the fiber has seven dimensions (two SU(2) and a mutual phase) but for the SU(3) 
subgroup of SU(4),both the base and manifold again reduce to four dimensions each. 
With z a pair of complex numbers, the non-trivial part of geometrizing SU(3) is thereby 
reduced to describing this four- dimensional manifold. Exploring this for the A^ = 3, 
n = 1 decomposition will be the content of the next section whereas the Appendix gives 
the alternative SU(4) rendering. 

3. Geometry of general SU(3) time evolution operator 

A general time-dependent three-level Hamiltonian may be written in terms of eight 
linearly independent operators of a three-level system. Such a Hamiltonian can also 
be written in terms of a subgroup of 15 operators of a four-level system. Before the 
time evolution operator is presented in the SU(3) basis in terms of a A^ = 3, n = 1 
decomposition, we will note that it can be rendered in a few alternative ways. 

First, a general time- dependent four-level Hamiltonian may be written as H{t) = 
X^jCjOj. Here Cj are time-dependent and Oj are the unit matrix and 15 linearly 
independent operators of a 4-level system that may be chosen in a variety of matrix 
representations. One choice used in particle physics are the so called Greiner matrices 
[361 [2n [22I [231 [2i] . Another choice consists of using ct, r, <t ® r and the 4x4 unit 
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matrix. Such a choice was discussed in [23 [2S] and will be used throughout this paper. 
As it stands, the above Hamiltonian describes a general four-level atom with 4 energies 
and 6 complex couplings. Note that only the three differences in energies are important. 
Restricting the 15 coefficients q to a smaller number allows this Hamiltonian to describe 
various physical Hamiltonians, forming different subalgebras of the su(4) algebra |25j . 
For example, if two of the six complex couplings are zero (levels 1 and 4 and levels 2 
and 3 of a four- level atom not coupled), then the Hamiltonian may be recast such that 
the operators involved belong to an so(5) subalgebra [25]. On the other hand, if levels 
2 and 3 are degenerate and level 4 is uncoupled from the rest, then the problem may be 
recast in terms of only eight operators belonging to the su(3) subalgebra of su(4). This 
is illustrated in Fig. [2] and is one of the systems of interest in this paper. 




Figure 2. Levels |2) and |3) couple equally to |1) and to |4), which are themselves 
coupled. The three complex coupling matrix elements and two energy positions define 
such an SU(3) system. 

Alternatively, after one arrives at the linear equation for the = 4, n = 2 
decomposition, one can represent the resulting vector in terms of six homogeneous 
coordinates. This is the so-called "Pliicker coordinate" representation for the SU(3) 
Hamiltonian. These coordinates as well as the alternative derivation are presented in 
the appendix. The N = 3, n = 1 decomposition will be the content of the rest of this 
section. 

Consider the Hamiltonian in the basis of the Gell-Mann lambda matrices [37] 
H{t) = J2i The A^ = 3, n = 1 decomposition consists of writing the time evolution 
operator in terms of a product of two matrices U = U1U2 where Ui is composed of a 
(2xl)-dimensional z, as explained in Sec. II. The equation that governs the evolution 
of z, equation ([6]), can be written in this 

= -iVf, - iF^^z^ + iV*z^Zf,; = 1, 2. (13) 

Here, the symbols used in defining z are defined as V = (04 — ia^, ag — ia-i), and 

_ ( ^3 + V^as ai - \ 
\ ai + ia2 —as + v^ag / 



Block sphere like construction of SU(3) Hamiltonians using Unitary Integration 8 

Using the transformation equations 2 = — 2:1,2 (-De**^)"^, = (De*'^)"^ and + 
|m2p + |m3p = 1 leads to the evolution equation for m = {rriir, rn2r, m^r, rriu, m2i, m^i)^: 
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rh, (14) 



which describes the rotation of a unit vector in a six dimensional space of the real 
and imaginary parts of fh defined by = m^r + ^^^i- In the above equations, 
D = {1 + l^ip + 1-22^)^^^ and i(j) = i{V*Zi, + Vuzl). The phase is real and determined 
only up to a constant factor. Since the real and imaginary parts of mg are not 
independently defined, the geometrical description of the base manifold for the = 3, 
n = 1 decomposition may be thought of as a point on the surface of a constrained 
six-dimensional unit sphere. 

The two constraints, namely |mip-|- |m2p + |m3p = 1 and the "phase arbitrariness" 
of 0, reduce the 6-dimensional manifold of the three-dimensional complex vector fh 
to a four- dimensional manifold in agreement with there being only four independent 
parameters in z.The first condition defines the base as a vector on an sphere while 
the phase arbitrariness serves as an additional constraint. The fiber, on the other hand, 
is an SU(2) block, evolving as a vector on S'^ Poincare-like sphere with a phase at each 
point, and a U(l) block that amounts to an extra phase. This is presented schematically 
in Fig.Q, as the product of three matrices of the evolution operator. 



1 


5^ 1 

^, 


F 
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'S'XU(I) ; 
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1 
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1 
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0' 





Figure 3. The base and fiber for tlie SU(3) group. Tlic first two factors give tlie 
base manifold, an spliere witfi a phase arbitrariness defined in the text. The fiber, 
described by the third matrix, is composed of a Bloch sphere and a phase associated 
with each of its points, and the second an extra phase represented by a vertical line. 

The alternative = 4, n = 2 decomposition in the appendix yields the equation of 
motion for = —z^/De^'^ in equation ( lA.llI) . Following equation ( ITTi) and equation 
( IT2|) . we see that for this case, the two remaining blocks of the time evolution operator, 
namely t/*^^"^-* and U^'^\ can be transformed into unitary matrices for SU(2). The fiber 
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evolves as vectors on two identical 5*^ Bloch-like spheres with a mutual phase, whose 
evolution is coupled to the base that evolves as a vector on an 5*^ sphere. This is 
illustrated in Fig. (jl]). Either decomposition can be used to study various physical 




Figure 4. The base and fiber for tlie SU(3) group via the iV = 4, n = 2 decomposition. 
The base again is given by an sphere as in Fig. ([3]). The fiber is composed of two 
identical SU(2) Bloch spheres plus phase, and an extra mutual phase between them. 
The four parameters each of base and fiber again account for all eight parameters of 
the SU(3). 



processes as will be discussed in the next section. 



4. Applications 

It is often desirable to control the time evolution of quantum states to manipulate 
an input state into a desirable output state. In [38l [39], the authors considered a 
Hamiltonian of the form Hq — fJ'S{t), where Hq is a free-field Hamiltonian and ^£{t) 
is a control field. To illustrate the "Hamiltonian encoding" scheme to control quantum 
systems, the authors considered a three-level system and studied stimulated Raman 
adiabatic passage (STIRAP), an atomic coherence effect that employs interference 
between quantum states to transfer population completely from a given initial state to a 
specific final state. This is done through a "counterintuitive" pulse sequence. Consider 
the Hamiltonian 

/ Giit) \ 
H(t)= Gi{t) 2A G2{t) . (15) 
V G2{t) ] 

Here G'i,2(i) = 2.5exp[— (t — ^1,2)^/^^] and A = 0.1. The initial population is in the 
upper state. For ti = r, t2 = and r = 3, it is seen that the two empty states are 
coupled first via 6*2 (^) and then the levels |1) and |2) are coupled through Gi. The 
dynamics of the populations reveal complete population transfer. A complete solution 
as per Section [3] was constructed for this model and the results are presented in Fig. [5] 
in total agreement with the results of [38] • 
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Figure 5. Population Pij = plotted as a function of time. The initial 

population in state |1) is completely transferred to |3). Both the unitary integration 
solution and the direct numerical solution [38] are plotted and they coincide at all 
times. 




Quantum control can also be achieved by understanding the nature of tunneling. 
The famous Landau-Zener formula O El [7] predicts the transition probability of the 
ground state of a two-level system when the energy levels adiabatically undergo a 
crossing. The study of level crossings has since been extended to multi-level systems. 
For example, in [ID], the authors considered a three- level atom to study population 
trapping by manipulating the phase acquired as a three-level system evolves under the 
influence of frequency modulated fields [H]. Such a frequency modulated field is given 
by 

E{t) = Eie-^['^i*+'^iWl + E2e-^['^2*+'^2Wl + c.c. (16) 
<^i{t) = Mi sin flit. (17) 

Here, c.c. stands for complex conjugation. The phase (fi{t) in the exponent can be 
written in terms of Bessel functions as [12] 

oo 

gAf,sinC,t_ ^ Jfc(Mj)e*^^^*. (18) 
fc=— oo 

For large values of Qj, the leading contribution for slow time scales would come from 
Jo{Mj). Hence, for large flj, the interaction Hamiltonian can be written as 

H,„i(t) = -d.(Ei Jo(Mi) + E2MM2). (19) 

Hence, for values of Mi^2 that are zeros of the zeroth-order Bessel functions, the 
interaction Hamiltonian is zero and population trapping is observed. Under this 
assumption, consider the full Hamiltonian under the rotating wave approximation, 

/ E,{t) Gi(t) \ 
H(t)= Glit) G2{t) . 

[ cm E,{t) ) 

Here, Ei{t) = Ai - MiQi cos(fiit + 6) and E^lt) = -A2 + M2Q2 cos(fi2t). Results are 
presented in Fig. [6], and for the parameter values ^1,2 = 1, Ai = — A2 = 10, 6* = and 
Gi,2 = 6, demonstrate the phenomenon of population localization discussed in [iO] . 
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Figure 6. (a) For Mi^2 = 7 and the other parameter values given in the text, there is no 
population trapping observed, (b) The energy landscape for Mi 2 — 30.6346 showing 
energy level crossing, (c) Population trapping is observed with Afi,2 = 30.6346 which 
corresponds to the tenth zero of the zeroth-order Bessel function. Note that the thick 
line is Pu and the thin line corresponds to P12. The results agree completely with 



As a final illustration of the unitary integration technique applied to three-level 
systems, let us consider the example discussed in |13]. Here, a three- level system is 
subject to strong fields and the correlation between the scattered light spectrum and 
the atom dynamics is discussed. The authors consider the Hamiltonian 



H(t) 



Gi{t) \ 
G2{t) 




(20) 






\G\{t) cm 

Here, Gi,2(^) = — Vi,2e~*''*. The time evolution of the states calculated as per our 
procedure in Section [3] is plotted in Fig. [7| for different values of the parameters. All of 
these results agree with those given in [l3j. Further features of the base and fiber will 
be presented at the end of the next section. 
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Figure 7. (a) Populations Pij = for S = 5, Vi = 2 and V2 = 1. Pii is given 

by the solid line and P12 is given by the thin line, (b) Same as (a), for S = 12. Note 
that Pi3 oscillates close to zero at all times, (c) Pij for 6 — 12, Vi — 1 and V2 — 2. 



5. Geometric phase for SU(3) group 

Many physical systems give rise to a measurable phase that does not depend directly 
on the dynamical equations that govern the evolution of the system, but depends only 
on the geometry of the path traversed by vectors characterizing the state of the system. 
This geometric phase is denoted by ■jg and is given by the integral [H], 

7, = 1 t/R .(n(R(t))|«VR|n(R(t))), (21) 

where the state evolution is governed by a set of internal coordinates that parameterize 
the Hamiltonian R(t), and Vr is the gradient in the space of these internal coordinates. 
This phase has been generalized to non-cyclic non-adiabatic evolution of quantum 
systems [15| B6| HTj. The purpose of this section is to present this phase in terms 
of coordinates on the Bloch sphere for two-level systems and extend it to three-level 
systems. 

In two-level systems, the time evolution operator is described by three parameters 
as described in Section [H Two of these parameters describe a point on the Bloch 
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sphere. Traversing closed loops on this Bloch sphere returns the quantum system to its 
initial state as described by the two parameters on the Bloch sphere but not the third 
parameter of an overall phase. Hence, general closed loops on the Bloch sphere do not 
correspond to closed loops in the space of the full unitary operator. This discrepancy in 
the phase between the initial and final state corresponds to the geometric phase given 
above and amounts to changes along the fiber at each point on the sphere. To formalize 
this, consider Ui, given by equation ([3]) as unitarized through the matrix b in Section [2l 
which for = 2, n = 1 takes the form 



By identifying cos | = (1 + |zp)~^ and sin |e~*'' = — z(l + |zp)^5j we get the usual 
description of the base manifold in terms of the angles < 6 < n and < e < 27r that 
are associated with the Bloch sphere, namely. 

In terms of the parameters 6 and e, the Hamiltonian H{t) = —a. a is given by 

H{t)=( -^i-^f^y (24) 

^ ' V - sin 9e'' cos 9 J ^ ^ 

equation governing the evolution of the fiber U2 has two terms. The first term is 
evaluated as 

UlH{t)U,= (^-^^ (25) 

which corresponds to the eigenvalues of the Hamiltonian. To evaluate the second term, 
consider the case whereby the vector on the Bloch sphere traverses a closed path defined 
by a constant 6. The second term is then given by 

rrl dU, _ ( -sin^f -|sin^e-M 

Integrating e from to 2tt yields 

f'\.rri dU, /7r(l-cosg) \ 

h ^d{-it) \ -7r(l -COS0) ) ' 

which is the correct formula for the geometric phase of a two-level system |44j . 

To extend this analysis to three-level systems, we consider the = 3, n = 1 
decomposition. The matrix Ui = Ui.b is now given by 



/(2) 



rZZ 




U, = I ^(^+^) ^ I , (27) 



where z is a complex column vector {zi, 2:2)^ and D = yl -\- |zp. Care has to be taken 
in assigning angles to elements of this matrix such that the transformation satisfies 
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two conditions: the Ui matrix should not depend on and the transformation must 
be commensurate with the definition of fh. To this effect, we transform z into polar 
coordinates: Zi = — tan y cos %-e*'^% Z2 = — tan ^ sin ^e**^^. These transformation 
equations imply that D = 



1+ z 



sin ^ sin ^e^^"^- 



mi 



( 



\ 



1 

sm 



and m3 = cos ye The Ui matrix is given by 
2^ -sin^ ^sin^2e*('^"'2) 

i{e-l-f-2) 



sin Y cos yc**^^^" 



1112 



- 2 sin^ ^ cos^ f 
2 ^ sin e-).e 



sin f cos f e-"i 



sin ^ sin 6'2e* 
1-2 sin^ ^ sin^ f 



sin Y sin ye" 



sin Y cos ye' 
sin Y sin yC*' 



cosf 



(2^ 



In the above equation, the range on the angles < < vr and < e-j < 27r 
are chosen so that the absolute value of each element of the time-evolution operator is 
positive [IE]. Hence Ui can be represented as two vectors on a sphere, at angles (6'i, ei) 
and (6*2, £2) respectively. This is represented in Fig. ([8]). Since the columns of a unitary 




Figure 8. The base manifold Ui is characterized by two sets of angles < 6',; < tt, 
< e.; < 27r which can be represented as two vectors with polar angles (6*1 , ei ) and 

(02, £2). 



operator correspond to normalized eigenvectors, we can consider the last column of 
the matrix above, = (— sin y cos ^6*"^, — sin y sin yc^^^, cos y)"^, and evaluate the 
so-called connection 1-form given by [19] 

A=-i{^\d\ij). (29) 

The Abelian geometric phase, given by 7^ = / ^ is evaluated to be 

If 9 

7g = -- y sin^ — {{dei + dt2) + cos02(c?ei - dt2)). (30) 

If the various angles are relabelled ei — *• —7 — a, €2 — ^ —7 + a, 61 ^ 29 and 6*2 2j3, 
the formula above agrees with [50] and [H] . The time-evolution operator above can now 
be used as in the case of SU(2) to evaluate the dynamic contribution / UlH{t)Ui and 
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the geometric contribution to the time evolution operator which is given by —i J UldUi, 
where dUi = ^dOi + ^de^ z = 1, 2. 

This description of the base manifold in terms of {6i, €{) can now be used to describe 
the dynamics of various physical processes. Fig. ([9]) represents the base manifold 
corresponding to the results in Fig. ([71). (6'i,ei) depend on all the parameters that 
define the system while (6*2,62) depend only on the ratio V1/V2. Also note that the 
maximum value of €2, corresponding to the maximum latitude traversed by the black 
curve, is inversely proportional to 6. Such observations can be used to control the 
dynamics of this system. 




Figure 9. The base manifold corresponding to tlie results in Fig. ((7|) for the three- 
level system of [43^. For the first column, Vi = 1, V2 = 2. The second column 
corresponds to T^i = 2, V2 = 2 and the third to Vi = 2, V2 = 1- The rows correspond 
to 5 = 1, S — 5 and 6 — 50. The thin black curve describes {9i, ei) and the thick red 
curve the set (6*2, £2)- 

6. Conclusions 

The ability to decouple the time dependence of operator equations from the non- 
commuting nature of the operators is the central feature of unitary integration and 
also characterizes the Bloch sphere representation for the evolution of a single spin. By 
doing so, the quantum mechanical evolution is rendered a "classical" picture of a rotating 
unit vector. For a two-level atom, the Bloch sphere representation along with a phase 
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completely determines the time evolution operator. In this paper, we have extended 
this program to deal with the time evolution operator belonging to the SU(3) group. 
This complements the work in [3l] for SU(4) Hamiltonians of two qubit systems. We 
have also extended the analysis of geometric phase to three-level systems by providing 
an explicit coordinate representation for the SU(3) time evolution operator. 



Appendix A. Alternative derivations for a general SU(3) Hamiltonian. 

Consider a three-level Hamiltonian written in terms of the Gell-Mann matrices [37] as 
H(t) = J^i^iaiXi. To exploit the fact that this Hamiltonian is a subgroup of four-level 
problems, it is represented in terms of the O matrices [25] as 

2-^02 + (ag - ^)03 + (2a3 + 2-^)04 + 0405 + a^Oe + ^a^O-j + 2a508 + 

aiOg + asOio + 2aiOii + 2a20i2 + 2a60i3 + 2a60i4 - 2ajOi^ + 2a^Oi&. (A.l) 

This embeds the Hamiltonian H{t) = J^i^iXi as a 4x4 matrix with zeros along the 
last row and column. In such a representation, the various entries of the Hamiltonian 
equation ([2]) are given by 

H'^^'^^ = ^agl^^^ + aiCTi + 020-2 + 030-3, (A.2) 
v3 

V = ^(04 - ^5)1(2) + l(a6 - ^7)0-1 (A.4) 

Writing z in the standard Clifford basis as z = ^z^I^'^^ — ^Y^i^iCTi, it follows from 
equation ([6]) that zi = iz2 and z^ = iz4^ and the equation reduces precisely to equation 
(IT^ . The geometry described in Section [3] can thus be derived from either of these 
decompositions of the time evolution operator. 

The SU(3) subgroup in equation (lA.ip is one among many SU(3) subgroups 
embedded in SU(4). Another choice corresponds to the Dzyaloshinskii-Moriya 
interaction Hamiltonian [121 113] and is also of interest because the 4x4 matrices now 
do not have a trivial row and column of zeros. In the two-spin basis, this Hamiltonian 
is given by 

Hit) = qO, = ai(02 + O3) + 2a2(Oi5 + Oie) + 2a3(Oi4 - O13) + 204(07 + On) 

i 

+05(05 + Oio) + 05(05 + O9) + 207(03 + O12) + ^(204 - Oi3 - O14). (A.5) 

The correspondence between the coefficients in terms of O and in terms of the A matrices 
is : ci = 0, C2 = oi, C3 = oi, C4 = 4o8/ a/3, C5 = Og, Cq = 05, C7 = 204, Cg = 207, Cg = oe, 
Cio = 05, cii = 204, C12 = 207, Ci3 = -203 - 2as/V5, Ci4 = 203 - 2o8/\/3, C15 = 2o2 
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and Ci6 = 2a2. Relabeling of the states 1— s>2, 2— s>3, 3-h>4 and 4 — > 1 expresses the 
Hamiltonian as 

-ff^"^"^^ = ^agl^^^ - a^cTi - a2(T2 - aia^, (A. 6) 

v3 

H^'^ = -^«8l(^) - ^asa„ (A.7) 
V = i(a6 - ^7)1(2) + l(a6 - ^7)0-1 (A.8) 

-^(os + «a4)cr2 - -(04 - ^asjcTs- 

If z is written in terms of the standard Clifford basis (I, —ia) as z = ^Z4l^'^^ — | I]^=i -^^i'^'i 
, it follows from equation ([6]) that zi = iz^ and ^2 = iz^- This is consistent with the 
parameter count that since the inhomogeneity V has only two free complex parameters 
(namely Vi = — ia-j and V2 = (24 — ia^)-, the complex z matrix should be composed 
only of two independent complex parameters, Zi and Z2. With the above analysis, 
equation ([6]) becomes for the pair of complex numbers 

]^z^ = ^Xf, - iF^^Zy + 2GyZuZ^] fi,u = 1, 2. (A. 9) 

Here X = {Vi/2, -iV2/2), G = {2V*,2iV2*) and 

■p _ ( ^^3 — V^^Cts 0-1 + ^0-2 \ 

\ —ai + ia2 —ia-i — \/^ia^ ) 
Paralleling the technique employed to solve an SO (5) Hamiltonian in [311 132] , 
we transform z into a complex vector m: = and ma = ^ such that 

|mip + |m2p + |m3p = l,with D = {1 + ^{\zi\^ + |z2n)"'"^^- This leads to the new 
set of evolution equations 

ias — ^f^ia^ oi + ia2 — og + ^a^ \ 
—ai + ia2 —ias — y/Sias as + ia4 m. (A. 10) 

ag + ia^ —as + ^04 J 

This can be written as an equation describing the rotation of the real and imaginary 

components of the vector fh = (mi^, m2r, m^r, f^u, fn2i, fT^u)^ 
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rh. (A.ll) 



Here, the coefficients Cj are written in terms of the coefficients a^, whose correspondence 
was given earlier in this section. Also note that = m^^ + D = (1 + |-2;2p)^ 
and = {V*Zi, + V^z*). Simplifying this leads to the equation i(p = — 2(X^z* — X*^z^) 
for the evolution of which is clearly real but determined only to within a constant. A 
little algebra yields for the effective Hamiltonian given by equation (ITT]) . 
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and for the effective Hamiltonian given by equation (1121) . the expression H^'^^ + (z^V + 
Vtz)/2. 

Another representation of the SU(3) subgroup of SU(4) Hamihonians is given by the 
so called "Pliicker coordinate" representation of the SU(4) group discussed in [311 [32] . 
For an arbitrary SU(4) matrix, the Pliicker coordinates are defined as a set of six 
parameters (P12, P13, P14, -P23, -^24, -P34) such that P12-P34 - -P13-P24 + -P14-P23 = and 



E \P^, 



1. They can be written in terms of the unit vector m and are given by 
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irriQ + ms 
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(A.12) 



The linear equation of motion for m translates into an evolution equation for P 
(P12, -Pi3, -P14, -P23, -P24, -P34) of the form iP = HpP. Here, Hp is given by 



H, 
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Vl 



PI 




(A. 13) 



where 
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In the above equation, Ojj± denotes Oi ± Oj. 
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